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ABSTRACT 

Lya photons that escape the interstellar medium of star-forming galaxies may be resonantly scat- 
tered by neutral hydrogen atoms in the circumgalactic and intergalactic media, thereby increasing 
the angular extent of the galaxy's Lya emission. We present predictions of this extended, low surface 
brightness Lya emission based on radiative transfer modeling in a cosmological reionization simula- 
tion. The extended emission can be detected from stacked narrowband images of Lya emitters (L AEs) 
or of Lyman break galaxies (LBGs). Its average surface brightness profile has a central cusp, then 
flattens to an approximate plateau beginning at an inner characteristic scale below ~0.2 Mpc (co- 
moving), then steepens again beyond an outer characteristic scale of ~ 1 Mpc. The inner scale marks 
the transition from scattered light of the central source to emission from clustered sources, while 
the outer scale marks the spatial extent of scattered emission from these clustered sources. Both 
scales tend to increase with halo mass, UV luminosity, and observed Lya luminosity. The extended 
emission predicted by our simulation is already within reach of deep narrowband photometry using 
large ground-based telescopes. Such observations would test radiative transfer models of emission 
from LAEs and LBGs, and they would open a new window on the circumgalactic environment of 
high-redshift star-forming galaxies. 

Subject headings: cosmology: observations — galaxies: halos — galaxies: high-redshift — galaxies: 
statistics — intergalactic medium — large-scale structure of universe — radiative 
transfer — scattering 



1. INTRODUCTION 

High-redshift star-forming galaxies are becoming an 
important probe of galaxy formation, reionization, and 
cosmology. Ionizing photons of young stars in star- 
forming galaxies ionize neutral hydrogen atoms in the 
interstellar medium (ISM), and each subsequent recom- 
bination has a probability of ~ 2/3 of ending up as a 
Lya photon (partridge fc Peebleslll967[) . After escaping 
the ISM, these Lya photons can be scattered by neutral 
hydrogen atoms in the circumgalactic and intergalactic 
media (IGM) , which tends to make the Lya emission ex- 
tended. In this paper, we present predictions for the ex- 
tended Lya emission associated with high-redshift star- 
forming galaxies from a radiative transfer model of Lya 
emissi on applied to a hy drodynamic cosmological simu- 
lation (jZheng et al.ll2010D . 

As a result of reprocessed ionizing photons, promi- 
nent Lya emission can be a characteristic of star-forming 
galaxies, which can be used to detect high-redshift 
galaxies. Galaxies detected through the strong Lya 
emission associated with them (e.g., from narrowband 
photometry) are dubbed Lyman-alpha emitters (LAEs) 
and an increasing number of such galax i es have been 
discovered (e.g., iHu & McMahonl 119961 iCowie & Hul 
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galaxies can also be detected from broadband photome- 
try with Lyman break techniqu e, which are termed Ly- 
man break galaxies (LBGs; e.g., St eidel et al.ll2003l) . As 
long as Lya photons reprocessed from ionizing photons 
can escape from the ISM, Lya emission is also expected 
to be associated with LBGs. Lya emission encodes useful 
information about star-forming galaxies a nd their envi- 
ronm ents, such as star formation rat e (e.g..lMadau et al.l 
119981 ). kinematics of ISM gas (e.g., iSteidel et al.l 120101). 
and neutral fraction of IGM gas (e.g.. iDiikstra et all 
2007). Lya photons usually experience complex radia- 
tive transfer processing (resonant scatterings) in media 
with neutral hydrogen atoms, which complicates the in- 
t erpretation of obser ved Lya emission properties. 

iZheng et al.l (|2010l ) (hereafter Paper I) present a phys- 
ical model of Lya emission from LAEs by solving Lya 
radiative transfer in the circumgalactic and intergalactic 
media. The Lya radiative transfer calculation is per- 
formed in a cosmological volume (100ft.~ x Mpc on a side) 
from a state-of-t he-art radiation- h ydrodynamic reioniza- 
tion simulation ()Trac et al.ll2008D . The calculation uses 
a 768 3 grid to sample the density, velocity, and tem- 
perature of the neutral hydrogen gas in the simulation 
and is applied to all the sources residing in halos above 
5 x 10 9 h~ 1 M ( T ) . Resonant scatterings enable Lya pho- 
tons to probe the circumgalactic and intergalactic envi- 
ronments (density and velocity structures) around star- 
forming galaxies. This leads to a coupling between the 
observed Lya emission properties and the environments. 
This simple physical model is able to explain an array of 
observed properties of z ~5.7 LAEs the Suba.ru/XMM- 
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Fig. 1. — Stacked narrowband Lya image for sources residing in 
halos of W 11 h- 1 M Q at z ~5.7. A mean algorithm is used to stack 
~440 sources. 



Newton Deep Survey (SXDS; iQuchi et al.|[2008h . includ- 
ing Lya spectra, morphology, apparent Lya luminosity 
function (LF), shape of the ultraviolet (UV) LF, and the 
distribution of Lya equivalent width. The selection im- 
posed by the environment dependent radiative transfer 
also intro duces interesting n ew features in the clustering 
of LAEs (jZheng et alJl201lL hereafter Paper II). 

In the above radiative transfer model, while the num- 
ber of Lya photons is conserved after they escape the 
ISM, the scatterings in the circumgalactic and inter- 
galactic media cause the Lya emission to spread spa- 
tially. Therefore, one generic prediction of the model is 
an extended Lya emission halo around a star-forming 
galaxy. Observationally, only a fraction of Lya photons 
can be detected for an individual source, those included 
in the central part of the extended Lya emission with 
high enough surface brightness (tip of the iceberg). The 
outskirts of the Lya halo with low surface brightness is 
typically buried in the sky noise. 

In this paper, we show that it is possible to detect the 
bottom of the iceberg by stacking the narrowband images 
of a large number of sources to suppress the sky noise. 
We present the predictions of the extended Lya emission 
in the stacked image from our radiative transfer model 
and discuss what we can learn from it. Our radiative 
transfer modeling is performed for sources at z ~ 5.7. We 
study properties of Lya surface brightness profile from 
stacked images in Section 2. In Section 3, we discuss 
the observational prospects. We summarize our results 
in Section 4. 

Throughout the paper, we adopt a spatially flat ACDM 
cosmological model for our calculations, with a matter 
density parameter f2 m = 0.28 and a Hubble constant 
h = 0.70 in units of 100kms _1 Mpc _1 . Distances are ex- 
pressed in eomoving units unless mentioned clearly oth- 
erwise. 

2. EXTENDED Lya EMISSION AROUND STAR-FORMING 
GALAXIES 

The Lya radiative transfer calculation in Paper I is 
performed in a simulation box of 100/i _1 Mpc (eomoving) 
on a side with the neutral hydrogen density, temperature, 
and peculiar velocity sampled with a 768 3 grid. The cell 
size is therefore 130.2ft, _1 kpc (eomoving), about the virial 



diameter of a 1.8 x lO lo /i _1 M halo. The star-forming 
region in a high-redshift galaxy, where Lya photons in 
our st udy originate, is only a few kpc (eomoving) across 
(e.g., iTaniguchi et all l2009f ). So the initial Lya emis- 
sion is treated as a point source at the halo center for 
the purpose of the radiative transfer calculation. The 
size of the pixel for collecting Lya photons is 8 times 
finer, corresponding to 16.3ft. -1 kpc (eomoving) or 0.58". 
With the cell resolution, the gas distribution is uniform 
inside the virial radius of small halos (below a few times 
10 10 /i _1 M Q ), which would have some effect on the Lya 
surface brightness profile. However, since Lya photons 
are initially emitted from a point source and the dynam- 
ics of the infall region around halos plays a significant 
role in determining the distribution of Lya photons (Pa- 
per I), the resolution we use is sufficient for obtaining 
the generic features in the extended Lya emission. We 
present a resolution test in § 2.5, and the main results 
presented in this paper would remain valid with improved 
resolution. Our calculation does not address the effects 
of galactic wind and dust, and we discuss these model 
uncertainties in 8 2.4 and S 3. 



2.1. Stacked Lya Image and Lya Emission at Fixed 
Halo Mass 

We start from the stacked Lya image for sources resid- 
ing in halos of fixed mass, 1Q H~ 1 Mq. In our model, the 
intrinsic Lya luminosity, which is the the total amount 
of Lya emission from the reprocessed ionizing photons, 
is tightly correlated with the halo mass, therefore this 
stacking is basically for sources at fixed intrinsic Lya lu- 
minosity or UV luminosity. 

Our radiative transfer model produces a three- 
dimensional array, recording the Lya spectra as a func- 
tion of spatial position on the sky. We construct the nar- 
rowband Lya image from this array, with a filter width 
similar to that in the z ~ 5.7 SXDS (see Paper I). Our 
simulation box has about the same area as SXDS with 
a redshift depth three times larger. So we are able to 
construct narrowband Lya images for three SXDS-likc 
fields. Each image corresponds to an ideal case with 
perfect continuum and sky subtraction. Examples of im- 
ages of individual sources can be found in Paper I (Fig. 4 
and Fig. 5). We stack all the source images together as 
would be done with the narrowband observation. For 
each source, the image includes the center and surround- 
ing pixels (up to a radius of ~ 10/i -1 Mpc). The center 
pixel is chosen to be the one that contains the halo cen- 
ter. In other words, the center corresponds to that in the 
UV band. It is evident that the surrounding pixels can 
include Lya photons from other neighboring sources. For 
all the sources in a narrow bin (0.16 dex) of halo mass 
around 10 11 /i _1 M Q (about 440 sources for each of the 
SXDS-like fields), we stack their images with their cen- 
ters aligned. 

Typically there are two algorithms to generate a 
stacked image, taking either the mean or the median. 
The stacked Lya image from the mean algorithm is 
shown in Figure [T] The mean algorithm has the ad- 
vantage of conserving the flux, even though it may be 
affected by outliers in the flux distribution, showing up 
as many small clumps in the stacked image. The me- 
dian stacked image has the advantage of insensitivity 
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to outliers. We find that at a fixed projected radius, 
the observed flux roughly follows a log-normal distribu- 
tion among different sources, which causes the surface 
brightness from the median stacked image to be lower in 
amplitude than that from the mean stacked image (see 
solid curves in Figure [2|). The stacking cases mentioned 
above correspond to an ideal situation of no photon noise. 
When total photon noise (source + sky) is added, the 
surface brightness profile from the mean stacked image 
remains the same and that from the median stacked im- 
age would change. In the limit of large noise (compared 
to the signal distribution), the flux distribution is com- 
pletely modified by the noise distribution, and the me- 
dian value is dr iven to be close to the mean (also see 
iPieri et all I2010D . In practice, the large noise limit is 
typical at large radii for the case considered here (see 
below) , so the surface brightness profile from the median 
stacked image approaches that from the mean stacked 
one. See the thin solid curve in Figure [5] In what fol- 
lows, we focus on presenting our results from the mean 
stacked image. We compute the surface brightness profile 
by averaging the flux in annuli at different radii, there- 
fore the effect of the grainy feature seen in the image is 
greatly reduced. 

The surface brightness profile from the stacked im- 
age shows a few interesting features. It has two char- 
acteristic scales that are associated with steep changes 
in the profile slope. The inner slope change occurs at 
i?i n ~ 0.1/i -1 Mpc (possibly upper limit, see § 12. 5j) and 
the outer one at i? ut ~ 1 Mpc. Inside the inner char- 
acteristic radius R ln , the profile appears as a central cusp, 
roughly following a steep power-law. Between the inner 
and outer characteristic radii, R{ n < R < i? ut, the sur- 
face brightness decreases slowly with increasing radius, 
close to a plateau. Beyond the outer characteristic radius 
i?out, the surface brightness profile steepens and drops, 
forming an extended tail. As shown later, the two scales 
have clear physical meanings. 

With the mean surface brightness profile, we compute 
the cumulative luminosity as a function of projected ra- 
dius (bottom panel of Figure [5]). It turns out that the 
luminosity from the stacked image does not converge at 
large radius. Because of the plateau in the surface bright- 
ness profile, the luminosity profile approximately follows 
R 2 , a trend similar to that from the global mean surface 
brightness (dot-dashed line). The global mean profile is 
computed by uniformly distributing the total Lya flux in 
the original narrowband image across the whole image. 
The plateau in the stacked surface brightness profile has 
a much higher amplitude than the global mean. As we 
show later, the higher signal is caused by the clustered 
neighboring sources. 

In the top panel of Figure [5J we also show the level 
of sky noise (dotted lines). The calculation is based 
on a 4-hour observation with the NB8I6 filter using 
the Subaru telescop^H, which is roughly t he setup for 
the s urvey of z ~ 5.7 LAEs in SXDS (jOuchi et al.1 
120081) . In the NB8I6 band, the sky (AB magnitude of 
20.4mag/arcsec 2 in a dark night) has a surface bright- 
ness of / s ky ~ 1-36 x I0 _15 ergs _1 cm~ 2 arcsec -2 . The 
mean sky photon count in a 2"diameter aperture from 

7 We use the Subaru Im aging Exposure Time Calculator at 
http: / /www. naoj.org/cgi-bin/img_etc.cgi 



a 4-hour exposure with the Subaru telescope reaches 
Asky ~ 7.1 x 10 6 . Therefore, the la sky noise in a 
2"diameter aperture corresponds to a surface brightness 
level of /sky/ ^/A s k y ~ 5.1 x 10 -19 ergs -1 cm -2 arcsec -2 , 
which is higher than the signal of interest on large scales 
(Fig. [2]). The noise per pixel is expected to be suppressed 
by viV with N sources stacked. The suppression works 
even better for detecting the signal at large radii, because 
the annuli used to compute the mean surface brightness 
profile can be much larger than the 2"aperture. Within 
each annulus, however, the noise in different pixels is cor- 
related, since the same pixel from the original image may 
fall into the annulus several times from centering differ- 
ent sources. The correlations need to be accounted for in 
estimating the noise level at each radius for the stacked 
profile. In practice, given the low surface brightness of 
the extended Lya emission and the high sky noise level, 
stacking the image before subtracting the continuum and 
sky may help to reduce errors caused by such subtraction 
and to increase the sensitivity. 

2.2. Decomposition of the Lya Surface Brightness 
Profile in the Stacked Image 

To understand the features in the stacked Lya surface 
brightness profile and the origin of the two characteristic 
scales, we perform a test by assigning artificial surface 
brightness profiles for individual sources. The image of 
each source (in halos above 5 x 10 9 /i -1 M@) is assumed 
to be a circular disk of radius R s with uniform surface 
brightness normalized such that the total luminosity is 
the intrinsic Lya luminosity from the source. That is, 
we replace the intrinsic point source with a uniform "see- 
ing" disk or (equivalently) modify it by a tophat point 
spread function (PSF). We then follow the same proce- 
dure as above to form the stacked image for sources in 
10 11 h~ 1 MQ halos and derive the Lya surface brightness 
profile. No radiative transfer is performed in this artifi- 
cial model, and we use it to illustrate the effect of source 
extent combined with source clustering. 

Figure [3] shows the series of surface brightness profiles 
with different "seeing" disk radii R s . For a small "see- 
ing" disk radius (R s = 0.08/i -1 Mpc), sources are close to 
point-like. The central spike (inside R s ) is clearly seen 
in the profile from the stacked image. Outside of this ra- 
dius, we still see signals, which are obviously contributed 
from neighboring sources as a result of clustering. 

The surface brightness profile has contributions from 
the central sources and from the neighboring clustered 
sources, analogous to the case of weak gravitational lens- 
ing around galaxies. In the spirit of the halo model (e.g., 
Seliak 2007I lBerlind fc Weinberg! 12002. ICoorav fc Shethl 
2002), the profile can be decomposed into two com- 
ponents, that from the individual source (the one-halo 
term) and that caused by the clustering (the two-halo 
term). Formally, we can write the "two- halo" term of 
the stacked surface brightness profile S 2 h(R.) as 

S 2h (R|M) = J dM'^L J d 2 H' 

[1 + w(R'\M, M')] E lh (R - R'|M'), (1) 

where dN/dM' is the surface number density of halos, 
w is the 2D two-point correlation function (essentially 



4 



Zheng et al. 



6 (arcsec) 



10 -ia 



o 

CD 

o 
u 

CO 



10- 



f-, 

CD 



PQ 

co 10" 20 



1 


10 100 


J. 1 1 1 1 


1 MINI 1 1 1 1 1 II 1 L 

[sky noise] 




— Mean Stacked 
V — Median Stacked 

— Median Stacked (w/ noise) 


-.Global 


Mean _ _ J_" ._!*_.- 




[sky noise]/10 v^i" 


"M h = 
i i 


lC'h-'Mg 

1 1 



0.01 0.1 1 10 

R (comoving h _1 Mpc) 



00 

Sh 

CD 



10 45 



10 4 



PS 10 43 

v 



10 



42 



1 


1 1 1 II II l i i i 1 1 1 i 


1 y I'l I it 




Mean Stacked 


■/ — 




Median Stacked / 






Median Stacked (w/ noise) / j 






// 






/// 

L //' 

intrinsic / // 










/// 






/// 






//£ 














E O /{■> 




i 


i i i I'm il 1 l 


- 

i 1 1 



0.01 0.1 1 10 

R (comoving h _1 Mpc) 



Fig. 2. — Surface brightness profile of extended Lya emission 
from stacked narrowband image for sources residing in halos of 
10 11 h~ 1 Mq. Top panel: the two thick solid curves are from mean 
and median stacked noise-free images, respectively. The associ- 
ated dotted curves indicate the scatter around the mean surface 
brightness, based on stacked images from three SXDS-like fields in 
our model. The thin solid curve is from the median stacked noise- 
added images, where the Poisson photon noise (source + sky) is 
computed for a 4-hour exposure with the Subaru telescope. The 
upper horizontal dotted line is the surface brightness level corre- 
sponding to the lcr sky noise in a 2" diameter aperture, assuming 
a 4-hour exposure with the Subaru telescope. It drops to the lower 
horizontal dotted line if 100 sources are stacked. The horizontal 
dot-dashed line is the global mean Lya surface brightness if the 
Lya flux from all the sources in each field is uniformly distributed 
across the field. Bottom panel: The cumulative Lye? luminosity 
distributions from the stacked images. The dashed line denotes 
the intrinsic Lya luminosity of one source residing in a halo of 
10 11 h^ 1 Mq. The dot-dashed line is the profile from the global 
mean surface brightness profile, which is oc R 2 . 

the angular correlation function), and Eih is the sur- 
face brightness profile of individual sources (the one-halo 
term). For simplicity, equation ([1]) assumes that de- 
pends only on halo mass, while in reality it depends on 
the environment around halos as well (see below). We 
emphasize that the one-halo term in the stacked profile 
is defined as the Lya emission from the central source 
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Fig. 3. — Lya surface brightness profile from (mean) stacked 
images with different assumptions about the extent of Lya emis- 
sion. Sources considered here reside in halos of 10 11 /! - 1 Mq. The 
thick solid curve is the prediction from the radiative transfer model, 
which is the same as that in Fig. [2] Other curves are for cases with 
no radiative transfer applied but with the Lya emission from a 
source artificially set to spread into a uniform "seeing" disk of ra- 
dius R s , which are to illustrate the effect of source extent combined 
with source clustering. The signal outside of this radius comes 
from the neighboring sources as a consequence of clustering. See 
the text. 

alone (i.e., contributed from all the sources to be center- 
aligned and stacked). Even though neighboring sources 
can have their Lya photons scattered into the one-halo 
scales of the central source, these photons still belong to 
the contribution from the two-halo term. 

In the limit of compact sources (small R s in the uni- 
form "seeing" disk case), £ih can be replaced with 
L(M')S D (K-K'), where L(M') is the intrinsic Lya lumi- 
nosity and 8d is the Dirac delta function. The two- halo 
term then becomes 



E 2 h(R|M) 



«M'^(M0[1 



w(R\M,M')}. (2) 



For the small "seeing" disk case in Figure [3] (R s = 
0.08/i _1 Mpc), the flattening below 0.3/i _1 Mpc is a sign 
of the halo exclusion effect, which leads to a zero signal 
in the 3D two-point correlation function and a constant 
in the angular two-point correlation function u>[3 

According to Equation ((2|), we expect to see a flat- 
tened profile on large scales, where id < 1. On these 
scales, which are out of the range of the plot, it is hard in 
practice to disentangle such a uniform background from 
the sky background. We therefore limit ourselves to the 
regime where w is still large. 

8 The bump around 0.3/i _1 Mpc can be understood as follows. 



R 2 



Roughly speaking, w(R) oc tl(R, Z)dZ with Zq 

the regime that R < r c , where r c ~ 0.2— 0.3h~ ^^Mpc is the halo 
exclusion scale and § is the 3D two-point correlation function (zero 
on scales below r c ). As R increases, the competition between a 
decreasing £ on average and an increasing path length of integration 
(i.e., a decreasing lower limit Zq) leads to a small bump in w around 
0.3h~ Mpc. If the seeing disk is comparable or larger than the 
size of the halo (e,g., the R 3 = 0.2h _1 Mpc case in Figure [3] for 
10 h~ x Mq halos), the bump is smeared out. 
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Fig. 4. — Mean Lya surface brightness profile from an individual 
source. The host halo mass for this source is 10 h~ 1 Mq. The 
thick and thin curves are the same source viewed along different 
directions. 

By varying the radius R s of the uniform "seeing" 
disk, we show the effect of the spatial extent of the 
Lya emission on the surface brightness profile from the 
stacked image (Figure [3]). As the "seeing" disk size in- 
creases, the two-halo term of the profile is smoothed 
on scales below the "seeing" disk size, which results in 
a decrease in the amplitude. In our setup, the uni- 
form "seeing" disk has a surface brightness of 1.4 x 
10" 19 (i? s /0.4/i" 1 Mpc)" 2 ergs" 1 cm" 2 arcsec" 2 for R < 
R s , which is by definition the one-halo term contribution. 
However, the small-scale surface brightness (in the one- 
halo regime) in Figure[3]decreases much slower than i?J 2 , 
as a consequence of the contribution from the smoothed 
two-halo term. For R s = 1.6/i _1 Mpc, the signal is al- 
ready dominated by the two-halo term on all scales. The 
approximate plateau seen in the stacked profile from the 
radiative transfer model (thick solid curve) is close to 
the case with a "seeing" disk size of R s ~ 1.6/i" 1 Mpc. 
The test suggests that the plateau feature between the 
two characteristic scales mainly originates from spatially 
extended, scattered Lya emission from clustered sources 
around the central source. 

The uniform "seeing" disk only gives us a rough idea on 
the features in the extended Lya emission. To be more 
accurate, we need to examine the Lya surface brightness 
profile of individual sources. 

As an example, we show in Figure |4] the surface bright- 
ness profiles of one individual source in the simulation 
box from the radiative transfer model. The two curves 
correspond to the profiles viewed along different direc- 
tions, and the difference reflects the environment de- 
pendent Lya radiative transfer, as discussed in detail in 
Paper I. The surface brightness has a steep drop below 
0.3/i -1 Mpc. Then it levels off until a cutoff is reached 
around 2-3/i" 1 Mpc. The inner surface brightness profile 
of one source roughly goes as R~ 2 and the luminosity 
increases logarithmically with radius. Most of the lumi- 
nosity comes from near the cutoff radius at ~2/i _1 Mpc, 
which is the reason why the plateau up to such a scale 
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Fig. 5. — One-halo and two-halo contributions to the Lya surface 
brightness profile. The dotted curve is the average profile from 
individual sources (one-halo term), and the dashed curve is the 
profile from source clustering (two- halo term). The solid curve is 
the sum of the two contributions, which would be the one from the 
stacked image. 

in the total profile is produced (see below). 

As discussed in Paper I, the shape of the surface bright- 
ness profile largely depends on the velocity field. Below 
0.3/z -1 Mpc, the contribution to the surface brightness 
profile mainly comes from the Lya photons that are last 
scattered in the infall region around the halo. On larger 
scales, where the profile levels off, the last scattering of 
Lya photons occurs in the IGM that undergoes Hubble 
flow. The cutoff radius corresponds to those photons that 
are the bluest after escaping the infall region. They can 
travel far away from the source before being redshifted 
into the line center to be significantly scattered in the 
IGM. The exact shape of the surface brightness profile 
also depends on the initial frequency of Lya photons (see 

As emphasized in Paper I and Paper II, Lya radia- 
tive transfer depends on circumgalactic and intergalactic 
environments (density and velocity structures) . Even for 
sources in halos of the same mass, there is a broad distri- 
bution of environments. Therefore, the one-halo term of 
the surface brightness profile is not identical for sources 
in halos of fixed mass or for the same source viewed along 
different directions (Figure 0|. The one-halo term in 
Equation ([T]) should have dependence on environment 
in addition to halo mass. To obtain the mean one-halo 
term for the stacked image and decompose the surface 
brightness profile, we perform Lya radiative transfer cal- 
culations separately for each source with halo mass in the 
lO 1:L ft. _1 M0 mass bin. We then stack the individual one- 
halo terms to obtain the mean. 

The dotted curve in Figure[5]is the mean one-halo term 
for sources in 10 11 h~ 1 M Q halos. The shape is similar to 
that seen in Figure |U The dashed curve is the two- halo 
term, inferred from subtracting the one-halo term from 
the total surface brightness profile. The two-halo term 
can be expressed as a convolution (Equation [1]), with 
the convolution kernel being the one-halo term of the 
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Fig. 6. — Surface brightness profile of extended Lya emission 
from stacked narrowband image as a function of observed Lya lu- 
minosity. All sources reside in halos of 1O 11 /i _1 M . Thin (thick) 
curves are for the faint (bright) half of the sources based on the 
observed Lya luminosity, which is contributed by the central, high 
surface brightness region. Dotted curves are the corresponding 
one-halo profiles. See the text for details. 

clustered star-forming halos of all masses. The shape 
of the one-halo term at different halo mass is similar to 
that shown in Figure [5j In the two-halo term, a plateau 
reaches out to a radius of l-2/i~ 1 Mpc, which is around 
the cutoff scale of the mass-averaged one-halo term. This 
is not a coincidence. The scale where the plateau in 
the two-halo term ends marks the spatial extent of the 
extended Lya emission of the clustered sources, since 
on scales larger than this the smoothing effect from the 
source PSF is small. As the two-halo term dominates on 
large scales, the above scale is just the outer characteris- 
tic scale i?out seen in the total surface brightness profile. 
The inner characteristic scale i?; n marks the transition 
from the one-halo term dominated regime to the two- 
halo term dominated regime. The realistic decomposi- 
tion of the surface brightness profile confirms the results 
of the test in Figure [3] that the origin of the plateau be- 
tween the two characteristic scales in the total surface 
brightness profile is the extended Lya emission. 

The environment dependent Lya radiative transfer 
couples the observed Lya emission with the circumgalac- 
tic and intergalactic environments. Many consequences 
of such a coupling are studied in Paper I and Paper II. 
For the stacked surface brightness profile, the environ- 
ment dependence also appears, as detailed below. 

At fixed halo mass, which corresponds to fixed intrin- 
sic Lya luminosity or UV luminosity in our model, the 
observed Lya luminosity has a broad distribution, re- 
flecting the distribution of environments. In our model, 
for each source we link together the pixels ab ove a surface 
bright ness threshold similar to that used in lOuchi et al.l 
(2008) to define the observed emission, and if none of 
the pixels around the central source exceeds the thresh- 
old, we simply take the flux in the pixel at the source 
position as the observed emission (see Paper I for de- 
tails). The observed luminosity is only a fraction of the 



extended Lya emission, coming from the central, high 
surface brightness part of the source. The one-halo term 
of the surface brightness profile is expected to be closely 
related to the environment. 

We divide the sources into two samples of equal size, 
according to the observed Lya luminosity. In Figure |6j 
we show the stacked surface brightness profile for the 
bright half and faint half, respectively. The dotted curves 
show the corresponding one-halo terms. We note that the 
fluxes inferred from the one-halo terms of the two cases 
are not necessarily the same. The reason is that the scat- 
tered emission is not isotropic, being sensitive to the den- 
sity and velocity structures around the source. Overall, 
the sample of sources with lower observed Lya luminos- 
ity has a shallower stacked surface brightness profile at 
small radius and there is no clear scale of the transition 
from the inner cusp to the approximate plateau. The re- 
sult implies that the stacked Lya narrowband image for 
star-forming galaxies with weaker observed Lya emission 
would appear to be less compact. 

2.3. Dependence on UV and Lya Luminosity 

So far we have focused on the Lya surface brightness 
profile for sources in halos of lQ ll h~ l MQ, The profile is 
expected to be a function of the halo mass. We show the 
dependence on halo mass in Figure both for halo in 
narrow mass bins (left panel) and above different mass 
thresholds (right panel). Since in our model the UV lu- 
minosity is tightly correlated with halo mass, the plot 
also represents a sequence of Lya profiles as a function 
of UV luminosity. 

The amplitude of the stacked profile increases with 
halo mass. On small scales this increase is largely a 
reflection of the fact that halos of higher masses host 
sources of higher Lya luminosity, while on large scales 
there is an additional contribution from stronger source 
clustering around more massive halos. 

On small scales (inside the inner characteristic radius), 
where the one-halo profile makes a substantial contribu- 
tion, sources in higher mass halos show a steeper profile. 
The inner characteristic scale i?i n , where the profile starts 
to level off, increases with halo mass. It implies that 
sources in higher mass halos have more extended Lya 
emission. The outer characteristic scale -R ou t, where the 
profile starts to steepen again, also increases with halo 
mass. As discussed above, this scale is an indication of 
the spatial extent of the Lya emission of the clustered 
sources. Massive halos are typically located in dense en- 
vironments, therefore the clustered sources around them 
tend to be in massive halos as well. On scales larger than 
-R ou t, the extended tail of the profile for higher mass halos 
is steeper. This is also a consequence of the higher clus- 
tering amplitude, since the surface brightness profile is 
proportional to (1+w) as in Equation ([1]), and the unity 
factor tends to flatten the profile for weakly clustered 
sources. 

Overall, the stacked Lya surface brightness profile is 
steeper for sources in higher mass halos, and the two 
characteristic scales increase with halo mass. On small 
scales where the central cusp dominates, if we approxi- 
mate the profile by a power law, we find that the power- 
law index decreases with halo mass, ranging from —0.5 
to —1.4 (from —0.5 to —1.7) for the mass bin (threshold) 
samples shown in Figure [JJ 
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The halo mass (UV luminosity) sequence does not have 
a priori selection based on observed Lya emission, so it 
applies to the case of LBGs. 

In practice, the brightness profile that can be easily 
measured in a narrowband survey is the stacked image of 
the identified LAEs, classified according to their appar- 
ent (observed) Lya flux. Figure [8] shows the stacked Lya 
surface brightness profile for LAEs as a sequence of the 
apparent (observed) Lya luminosity, for both luminosity- 
bin samples (left panel) and luminosity-threshold sam- 
ples (right panel). 

The apparent Lya luminosity couples to the circum- 
galactic and intergalactic environment around sources as 
a result of the Lya radiative transfer, which imposes a 
strong selection effect for sources that are identified as 
LAEs. On small scales (inside Ri n ), the stacked profile 
for LAEs appears to be steeper than those seen in cases 
in Figure [JJ for sources without Lya selection. In Fig- 
ure [SJ the central cuspy profile has a slope of— 0.9 1.8 

(—1.1 - —1.8) for luminosity bin (threshold) samples, 
with steeper slopes for higher observed Lya luminosities. 
The amplitude increases with the apparent Lya luminos- 
ity. On scales above the inner characteristic scale, the 
amplitude also increases with Lya luminosity, but the 
dependence is much weaker. This weak dependence is 
a consequence of the weak dependence of the projected 
clustering on apparent Lya luminosity, as discussed in 
Paper II. There is also a trend that the two characteris- 
tic scales increase with the apparent Lya luminosity, but 
it is much weaker than that seen in Figure [7] for LBGs. 

2.4. Dependence on the Intrinsic Lya Line Profile 

In all our discussions above, the intrinsic Lya line pro- 
file is assumed to be Gaussian with the width deter- 
mined by velocity dispersion in the halo (Paper I) . More 
specifically, the standard deviation of the Gaussian pro- 
file is set to equal the ID thermal velocity dispersion for 
gas in the halo, a = Sl.9[M h /{10 w h- 1 M Q )] 1 / 3 kms~ 1 . 
As discussed in Paper I and Paper II, the intrinsic line 



profile is one of the main uncertainties in our radiative 
transfer modeling. The spatial extent of the Lya emis- 
sion depends on the initial line profile. For example, in 
the extreme case where all the Lya photons are shifted 
far to the red from the line center within the emitting 
galaxy (e. g., as a possible conseque nce of a strong galac- 
tic wind; Di ikstra fc Wvithel I2010D . they would not be 
further scattered in the surrounding medium and the 
Lya emission would not extend beyond the size of the 
galaxy. 

To test the effect of the initial line profile on the spatial 
distribution of Lya photons, we perform radiative trans- 
fer calculations with different initial values of Lya wave- 
length for an individual source in a halo of 10 11 ft, -1 Mq 
(the same source as in Figure 0|) . Figure |9] shows the 
Lya surface brightness and cumulative luminosity pro- 
files for cases with different initial line shifts (in units of 
the thermal velocity dispersion a mentioned above). 

On average, bluer photons travel farther in the IGM 
before they redshift into the line center to encounter sig- 
nificant scatterings. Therefore, the Lya emission would 
dilute to a larger spatial extent. This is clearly seen in 
Figure [9] - the spatial distribution becomes more and 
more extended as the initial Lya wavelength shift varies 
from — 3er to —9a. If we start with red Lya photons, the 
final Lya surface brightness profile becomes more com- 
pact. In Figure [9j the intrinsic Lya luminosity for all 
curves is the same. The surface brightness curves for 
initially red photons appear to be of lower amplitude, 
because they have a central spike (within 0.02/i _1 Mpc) 
that is not shown in the left panel. The contribution 
of the central spike to the total luminosity can be fig- 
ured out from the cumulative luminosity curves shown 
in the right panel. Overall, the Lya emission appears 
to be close to point-like if photons start from the far 
red side of the line center. These red photons can still 
have a probability of experiencing resonant scatterings 
in the infall region around the halo and being shifted to 



Zheng et al. 



CM 

I 

o 

0) 

w 

CJ 
U 

CO 



CSD 

u 

CQ 
m 



10 



6 (arcsec) 
1 10 100 



10 



10 



-17 



1 I I I III 1 — I I I I 1 1 1 1 I I I I III 1 — c 



-18 




logL, = 42 



apparent 

j I i I i r^i 1 1 1 1 



-19 

0.01 0.1 1 10 

R (comoving h _1 Mpc) 



10 



o 

W 

o 
u 

cd 



9 (arcsec) 
1 10 100 



-17 



10 



-IB 



S-i 



i mns — i — i 1 1 1 1 1 1 — i 1 1 1 1 1 r 

.logL „ nt >43.6 

{ to apparent 



DQ 

in 10 -19 

0.01 




logL appa . e „t>42.0 
i I i I i ublllu 



0.1 1 10 

R (comoving h _1 Mpc) 



Fig. 8. — Dependence of mean Lya surface brightness profile on the apparent (observed) Lya luminosity. Left: Profiles for sam- 
ples of galaxies in different bins of apparent Lya luminosity. From bottom to top curves, apparent Lya luminosity increases from 
log[L apparcn t/(erg s" 1 )] = 42.0 to 43.6, with an interval of 0.2 dex between two adjacent curves. Right: Similar to the left panel, but 
for samples of galaxies above different thresholds in apparent Lya luminosity. 



10" 



6 (arcsec) 

10 100 



17 



g io- 18 

o 

* 10-19 
w 

i 

| 10" 20 

w 10" 21 

<u 10~ 22 



OQ 

m 10 -23 

0.01 




M h =10 11 h- 1 M Q 



J I I I I LHj 



0.1 1 10 

R (comoving h _1 Mpc) 



6 (arcsec) 
1 10 100 



10 43 



ffl 10 42 



A 

| — ' 10 41 



i i 1 1 in 1 — i i i 1 1 1 n 1 — i i i 1 1 1 ii — i — r 




.•"init. A shift. 
..** -9a - 

-6a - 

+3ct = 

M-^IO 11 !!-^ --+§<' : 

h — — +9(7 - 

I I I I I T 



0.01 0.1 1 10 

R (comoving h _1 Mpc) 



Fig. 9. — Dependence of Lya emission distribution profile on initial wavelength of Lya photons. The left and right panels are for the 
surface brightness profile and cumulative luminosity profile, respectively. For each curve, the initial Lya photons are set to have a single 
wavelength, denoted in either panel as the difference from the line center in units of the fD thermal width of the halo gas, which is about 
a = 69kms _1 for Mh = lO 11 ^ - X M0. The initial Lya luminosity for all cases is normalized to be the same. 



the blue side. The surface brightness profile on scales 
< l/i _1 Mpc is expected to be contributed by these pho- 
tons. The probability of such scatterings depends on the 
magnitude of the shift and the infall velocity around the 
halo. For sources in 10 11 h~ 1 M ( ? ) halos, even with a +6a 
initial line shift (thin dashed curves in Figure EJ), the 
flattening part in the surface brightness profile can still 
be clearly seen (left panel), although the majority of the 
observed Lya photons are from the central region (right 
panel). We verify that for a +3.5(7 initial line shift (not 
shown in the plot) the contribution from the central re- 



gion to the total Lya luminosity is <40%, less than that 
from the extended emission. If we take this as a thresh- 
old for the profile to change from point-like to extended, 
it would correspond to ~ 250km s . 

If the wavelength shift of the initial Lya photons is in 
the range of ±3<7, the resultant surface brightness distri- 
bution is not sensitive to the initial Lya wavelength, as 
shown by the solid curves in Figure [9] This is easy to 
understand: photons starting at a wavelength shift ±3<r 
undergo many scatterings in the core, leading to a distri- 
bution that loses memory of the original wavelength and 
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to a nearly constant surface brightness profile. But pho- 
tons with a larger initial wavelength deviation are not 
scattered enough times to be redistributed in this way 
before they escape. Therefore, the predicted Lya profile 
should be reliable except when many photons are emit- 
ted from the central galaxy with a wavelength shift of 
several times the velocity dispersion. 

To eliminate most of the extended Lya emission 
around the galaxies, the initial line shifts (after photons 
escape the ISM) need to be larger than 250km s _1 for ha- 
los of 10 11 h~ 1 M(7 ) . Galactic winds can have the effect of 
shifting Lya emission toward re d (e.g.. iVerhamme et al.1 
120061: iDiikstra fe Wvithd I2010D From observations of 
z ~ 3 galaxies. iSteidel et al.l (|2010l ) find that the in- 
terstellar absorption line shifts with respect to the Ha- 
defined galaxy systemic redshift by 165 ± 140km s _1 
(mean ± standard deviation) for galaxies of baryon mass 
(stars plus cold gas) ranging from ~ 1 x 10 10 M Q to 
~ 3 x 10 11 M Q . The uncertainty of velocity measured 
for individual galaxies is estimated to be ~ 130km s _1 . 
If the wind is confined in a region not far from the ISM, 
the above observational results indicate that the initial 
shift in Lya emission could be above ~ 300km s" 1 . How- 
ever, even in such a case, the shift may not be suffi- 
cient for Lya photons to completely decouple from the 
circumgalactic and intergalactic environments, depend- 
ing on the spatial distribution of the wind. Observa- 
tionally, the outflowing neutral gas of the wind is not 
isotropic, usually collimated in a bipolar fashion (e.g. , 
Bland fc Tullvl 119881: IShonbell fc Bland-HawthorrilfoM 
Veilleux fc Rupka 12002). It is likely that Lya emission 
escaping in directions other than the wind would gain 
only a smaller redward shift, and the source would still 
become extended in the end. The wind effect on the spa- 
tial distribution of Lya emission is also expected to be 
less significant for sources before reionization than those 
after reionization, since the higher neutral density boosts 
the scattering optical depth. 

It is also possible that the wind could reach sc ales much 
larger than the ISM (e.g., ISteidel et all 120101 ) . In such 
a case, the initial line shifts (after photons escape the 
ISM) is not much affected by the wind, but the wind 
alters the distribution and dynamics of neutral gas in 
the circumgalactic region and therefore the details in the 
Lya radiative transfer process. While the surface bright- 
ness profile of Lya emission may change accordingly, the 
extended Lya emission should remain as a generic fea- 
ture. Such a possibility deserves further investigations 
using high resolution galaxy formation simulations with 
a wind prescription. 

A change in the initial Lya line profile can lead to a 
change in the observed Lya luminosity and therefore the 
observed Lya luminosity function. However, as discussed 
in Paper I, a combination of changes in the star forma- 
tion prescription used, the stellar IMF adopted, and the 
dust extinction is able to keep the predicted Lya lumi- 
nosity function to match the observed one. On the other 
hand, the change in the initial Lya line profile can al- 
ter the strength of the coupling between the observed 
Lya emission and the circumgalactic and intergalactic 
environments. This would show up in the clustering of 
LAEs (Paper II), which provides an interesting way to 
constrain the initial Lya line profile and the effect of 
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brightness profile. Shown here are average (one-halo term) profiles 
from stacking narrowband images of individual sources residing in 
halos of 10 1 h~ 1 Mq. The solid curve is obtained by perform- 
ing Lya radiative transfer calculation with gas properties (density, 
temperature, and peculiar velocity) represented with a 768 3 grid 
(default case in this paper), while the dashed curve shows the case 
with a lower resolution grid (384 3 ). For each case, the thin and 
thick arrows mark the half and full size of a grid cell, respectively. 
Note that the image pixels to collect Lya photons is much finer 
than the grid cell. Calculation with the higher resolution grid tends 
to predict more extended Lya emission. 

galactic wind. 

2.5. Effect of Grid Resolution 

For the main results in this paper, our Lya radiative 
transfer calculation is performed in a 100 3 /i~ 3 Mpc 3 sim- 
ulation box, where gas properties (density, temperature, 
and peculiar velocity) are sampled with a 768 3 grid. All 
Lya photons are originally emitted from point sources 
and the scatterings of Lya photons inside individual cells 
are also followed. As shown in Paper I, the dynamics of 
the infall region around halos plays a significant role in 
determining the distribution of Lya photons, so we ex- 
pect that the grid resolution in our calculation is able to 
predict the generic features in the extended Lya emis- 
sion. We perform a resolution test to show that this is 
the case. 

It would be ideal to do such a test with a higher res- 
olution grid. However, the available reionization simu- 
lation outputs prevent us from doing this. Therefore, 
we degrade the default 768 3 grid used in this paper to 
a lower resolution 384 3 grid by binning the gas density, 
temperature, and peculiar velocity. We then perform 
the Lya radiative transfer calculation with such a grid 
for 10 11 /i _1 M Q halos and compare the resultant average 
Lya surface brightness profile (the one-halo term) with 
that with the 768 3 grid. 

The comparison is shown in Figure fTOl For each resolu- 
tion case, the thin and thick arrows mark the half and full 
size of a grid cell, respectively. For a source located at the 
center of a cell, the initial Lya photons see a uniform gas 
distribution within a radius of about half cell size. How- 
ever, the final Lya surface brightness distribution inside 
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this resolution radius is determined by the gas properties 
not only in this central cell but also in the surrounding 
cells. The reason is that Lya photons propagating into 
surrounding cells have chances to be scattered back into 
the central cell. That is, the surface brightness profile 
is not a local quantity as a result of radiative transfer. 
Therefore, the small scale surface brightness profile (in 
particular the slope) is not completely a result of the grid 
resolution, as indicated by the fact that the powerlaw-like 
inner profile naturally extends to scales a few times larger 
than the resolution scale. The transition to a plateau is 
not an artifact of grid resolution, either. 

With the low resolution grid, the general features in the 
Lya surface brightness profile remain the same, a steep 
drop on small scales and an extended plateau on large 
scales before a sharp cutoff. However, the amplitude of 
the profile on small scales is about a factor of two higher 
with the low resolution grid, and it is compensated by 
a slightly lower amplitude on large scales. Therefore, a 
lower resolution has the effect of making the emission 
more concentrated. 

If the grid were a factor of two finer than the de- 
fault grid, we expect the predicted small scale surface 
brightness profile to be lower and the difference should 
be smaller than a factor of two. The Lya emission would 
become more extended. We therefore expect that the 
generic prediction of the extended Lya emission would 
be strengthened with an improved grid resolution. 

Even though for the one-halo term of the surface 
brightness profile the overall cusp-to-plateau behavior 
does not depend on the grid resolution, the transition 
scale decreases for a higher resolution. Figure [TU] seems 
to show that the transition happens at a scale about twice 
the cell size, indicating a resolution effect. Although cells 
around the source cell contribute to the determination of 
the profile of the central cusp, the cuspy profile itself 
is inevitably expected to also depend on the exact den- 
sity, temperature, and velocity distributions of neutral 
gas within the source cell. The effects discussed in § 12.41 
add further uncertainties in the cuspy profile. 

Despite the above uncertainties, our prediction on the 
existence of the extended Lya emission around high red- 
shift star-forming galaxies remains valid, as long as Lya 
photons encounter scatterings in the circum-galactic and 
inter-galactic media. The predicted two characteristic 
scales in the surface brightness profile from stacked im- 
ages should also remain valid. The outer scale, which 
is mainly determined by the one-halo plateau region, is 
less sensitive to the exact central cuspy profile. The inner 
transition scale, on the other hand, depends on the am- 
plitude and slope of the central cusp, which may not be 
properly resolved in our model. If the grid resolution is 
the only source of uncertainty, the inner transition scale 
of 0.1-0.2 Mpc (comoving) seen in the stacked profiles in 
§[2] should be regarded as upper limits. 

3. DISCUSSION 

The Lya emission from star-forming galaxies is ex- 
tended because of resonant scattering in the circum- 
galactic and intergalactic environments. In the stacked 
narrowband Lya image, the change from a declining to 
nearly flat surface brightness profile at a characteristic 
scale of a few tenths of Mpc (possibly an upper limit, see 
§ 12. 5p . and that from the flat profile to a declining pro- 



file at a few Mpc are evidence for the resonant scattering 
nature of the extended Lya emission. 

According to our radiative transfer modeling, the vis- 
ibility of Lya emission from star-forming galaxies de- 
pends on the circumgalactic and intergalactic environ- 
ments around sources, and only the central, high surface 
brightness part can be detected as LAEs. The extended 
Lya emission around the star-forming galaxies is pre- 
dicted to exist (as long as Lya photons escape from the 
ISM), no matter whether they can be detected as LAEs 
or not. This implies that the extended emission can be 
detected from both LBGs and LAEs. In practice, to de- 
tect the extended Lya emission around LBGs, we need to 
know their redshifts and take their Lya narrowband im- 
ages. Extended Lya emission from LAEs is more readily 
detected, since LAEs are typically discovered in narrow- 
band surveys. 

Our radiative transfer modeling is performed for z = 
5.7 star-forming galaxies. The generic picture should 
remain valid at other redshifts - as long as Lya pho- 
tons encounter significant scatterings in the circum- 
galactic/intergalactic media, the Lya emission would be- 
come extended. However, the physical conditions in and 
around star- forming galaxies evolve with redshift. For 
example, compared to z = 5.7 galaxies, we expect a 
lower neutral hydrogen density around z ~ 3 galaxies, 
as a result of the expansion of the universe and the evo- 
lution of the UV background. The Hubble expansion 
rate decreases, a change that may partially compensate 
the decrease in neutral density for Lya scattering. The 
extended Lya emission (e.g., the surface brightness pro- 
file) at z ~ 3 may not follow the details of our 2 = 5.7 
prediction, although it is expected to show up. We re- 
serve the investigation of the redshift dependence of the 
extended Lya emission for future work. We do not in- 
clude dust in our model. The effect of ISM dust on the 
Lya lin e depends on whether dust distribu tion is clumpy 
or not (|Neufeldlll99lH H ansen !TOh|[200ll . Optical, UV, 
and Lya observations of local star-forming galaxies show 
evidence that ISM kinematics and geometry play a more 
signifi cant role than dust in affecting the Lya emission 
(eV iGiavalisco et aLlfLoM lK^el[200a lAtek et al.l[2rM 
2009). We also expect the dust content to decrease to- 
wards high redshift (e.g., z ~ 6), and there is observa- 
tional evidence for low dust content in z > 6 galaxies 
(e.g.. lBouwens et al.ll20"l0t lOno et alJl2010f ). It is likely, 
however, that the existence of dust can lead to varia- 
tions in the amount of Lya photons escaping from the 
ISM, and thus add additional scatter in the extended 
Lya emission among individual galaxies. Keeping all the 
above caveats in mind, we discuss the current status of 
observing the extended Lya emission. 

Observations have long been suggesting the existence 
of extended Lya emission. Based on Hubble Space Tele- 
scope (HST) broad-band and New Technology T el escope 
(NNT) narrow- band images, IMoller fc WaTrenl (fl998h 
found that three LBG-like galaxies near a z = 2.8 quasar 
are more extended in Lya emission (half-light radius of 
0.5"-0.7") than in the continuum (half-light radius of 
0.13"). Observations with NNT and Very Large Tele- 
scope (VLT) of a larger sample of LAEs (~ 70) in quasar 
fields at z ~ 3 furth er strengthen the c a se for extended 
Lya emission (e.g.. iFvnbo et all 120011 120031 : also see 
lLaursen fc Sommer-Larsenl 120071 for comparing one of 
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their objects with the results from a model wi t h Lya 
radiative transfer). In particular, iFvnbo et al.l (|200lD 
found the size of Lya emission to be 0.65" (FWHM) in the 
stacked narrow-band images of seven LAEs. Although 
the spatial extent of the Lya emission from the above 
observations is small compared to what we predict, ow- 
ing to the surface brightness detection limit, it is clear 
that Lya emission is more extended than the underlying 
UV emission. 

There are efforts to use stacked narrowband images 
in deep s urveys to study t he co mpactness of the Lya 
emission. lHavashino et al.l (|2004f ) f ind that 24 spectro - 
scopically confirmed z ~ 3 LBGs in lSteidel et al.l (|2003l ) 
fall within the redshift range of their deep narrowband 
image observed with Subaru telescope (with only two of 
them being identified as LAEs). The UV luminosity of 
these LBGs ranges from -20 m ag to -22.5 mag (bas ed on 
the observed IZab rnag nitude; ISteidel et all 12003). Af- 
ter removing two gigantic Lya blobs, the composite nar- 
rowband (continuum subtracted) image of the 22 LBGs 
shows Lya emission that extends at least to 5" (about 
0.1/i -1 Mpc comoving; see their Fig. 8). The surface 
brightness drops from ~ 3 x 10 -18 ergs -1 cm -2 arcsec -2 
at ~ 1" to a few times ~ 3 x 10 -19 ergs -1 cm -2 arcsec -2 
at ~ 4", corresponding to a power-law index of ~ —1.7. 
They further note that 13 LBGs with no individual de- 
tection of extended Lya emission also show extended 
Ly a emission on their co mposite image. The findings 
of lHavashino et al.l (J2004) are in line with our predic- 
tion that extended Lya emission exists for LBGs, and 
the slope of the surface brightness profile is also in broad 
a greement with our pr ediction for UV bright sources. 

ISteidel et al.l ([201 ll ) present the result of stack- 
ing analysis of deep narrow-band observation of 92 
UV continuum-selected star- forming galaxies at 2.2 < 
z < 3.1. El Some sources are the same 
ones as in lHavashino et al.l ([2004) but with much 
deeper narrow-band observation. The stacked narrow- 
band image shows diffuse Lya emission around star- 
forming galaxies, with surface brightness drops from 
a few x 10 -18 ergs -1 cm -2 arcsec -2 at the center to ~ 
10 -19 ergs -1 cm -2 arcsec -2 at a radius of ~ 10" (~80 
proper kpc). The extended emission exists in the stacked 
image of sub-samples - LBGs, LAEs, and Lya absorbers 
all show qualitatively similar diffuse emission. The over- 
all surface brightness amplitude of the extended emis- 
sion increases with the central Lya emission. Although 
the redshift difference prevents us from a direct com- 
parison between the observation and our prediction, the 
detected diffuse Lya emission supports our generic pic- 
ture that resonant scatterings in the circum-galactic and 
intergalactic media causes Lya emission originated from 
star formation in the central galaxy to become extended. 
If we make a naive comparison, the Lya emission is de- 
tected up to the inner characteristic scale in the surface 
brightness profile (e. g., Figures [7| and Et , not reaching 
the plateau yet. In ISteidel et all (pfflll). the clustered 
sources are intentionally masked before image stacking, 
which may change the shape of plateau if it could be 

9 It is exciting that our generic prediction (if dated from the 
initial submission of our paper) started to be tested or verified in 
merely three months, although many details remain to be investi- 
gated. 



detected. We plan to extend our model to z ~ 3 star- 
forming galax ies and make direct comparison with the 
o bservation in Stcidcl et al.l (|2011[) . 

IQno et all (120101 ) use stacked multiband images of 
LAEs at z ^5.7 and 6.6 in the SXDS to investigate the 
stellar population of LAEs. For the narrowband stacked 
images, they find that the Lya fluxes level off for aper- 
tures larger than 5" in diameter. Such a behavior in the 
surface brightness profile is predicted in our model as the 
transition to the two-halo ter m dominated reg ime. It is 
likely that the feature seen bv lOno et al.l (|2010D is caused 
b y the clustering of Lya e mitting sources. 

Finkclst ein et al.l (|2011l ) detect Lya emission from 
three spectroscopically confirmed z — 4.4 LAEs in the 
HST/ACS F658N narrowband imaging data. They find 
that Lya emission appears to be more extended than 
the UV emission. However, the half-light radii of Lya 
emission are quite small, in the range of 0.1"-0.2". One 
thing to notice is that the observation is not deep enough 
to reach the sensitivity of detecting the much larger ex- 
tended Lya halo predicted by our model. This is clear 
from the fact that the surface brightness level in the 
HST narrowband observation turns out to be above a few 
times 10 -16 erg s -1 cm -2 arcsec -2 . It is also supported by 
the fact that the three objects have significantly higher 
Lya fluxes detected from ground-based narrowband im- 
age (with larger telesc ope apertures and redu ced sen- 
sitivity to read noise; IFinkelstein etaLl |20TTI) . With 
HST/WFPC2 nar rowband images of eight z = 3.1 LAEs, 
Bon d et al.l (|2010D find half-light radii of Lya emission to 
be in a similar range, 0.15"-0.30". The surface brightness 
level (a few times 10 -17 ergs -1 cm -2 arcsec -2 ) reached 
by the observation is still too high to probe the extended 
Lya emission predicted by our model. 

Extended Lya emission can also be detected through 
deep spectroscopic observations. IRauch et al.l (2008) 
conduct a long-slit search for 2.67 < z < 3.75 low sur- 
face brightness Lya emission from a 92-hour long ex- 
posure with the VLT, reaching a la surface brightness 
detection limit of 8 x 10 -20 erg s -1 cm -2 arcsec -2 in a 1 
arcsec 2 aperture. They find 27 objects with possible Lya 
emission lines. These o bjects are typically fainter than 
the LAEs and LBGs in lSteidel et al.l (|2011h . About one 
third of them show clear signatures of extended emission. 
If the surface brightness profile at large radii is approxi- 
mated by a power-law, the slope is generally in the range 
of —1 to —3 (their Figure 17). In the stacked surface 
brightness profile (their Figure 20), the emission clearly 
extends to > 4" (about comoving 90/i -1 kpc), dropping 
from ~ 8 x 10 - erg s -1 cm -2 arcsec -2 at a radius of l"to 
~ 1 x 10 -19 ergs -1 cm -2 arcsec -2 at large radii. The pro- 
file roughly follows a power law with an index in the range 
of —1.5 to —2.0. This is in the ballpark of the predicted 
extended emission, although a direct comparison is not 
possible because of the redshift difference. The spatially 
resolved spectra of the extended emission remain as an 
interesting aspect to be explored in our model. 

The extended Lya halo from Lya scatterings pre- 
dicted by our model shares some similarities with the 
one around a quasar before r eionization descr i bed i n 
ILoeb fc Rvbickl (fl999l) . Unlike lloebfc Rvbickil (fl99ll . 
who assume a uniform, zero temperature IGM under- 
going Hubble expansion, we use a realistic distribution 
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of gas density, temperature, and velocity around a star- 
forming halo from cosmological reionization simulation. 
We expect that extended Lya emission around star- 
forming galaxies also exists before reionization, which 
merits further study. However, the Lya surface bright- 
ness is expected to be proportional to the source lumi- 
nosity and to not depend on the neutral fraction in the 
IGM as long as the majority of photons are scattered, 
while it should decline with redshift as (1 + z)~ 4 , making 
the detection increasingly difficult at higher redshift. 

Our radiative transfer calculation is done for sources 
residing in halos above 5 x 10 9 /i _1 M Q . Star formation 
in lower mass halos, down to ~2 x 10 8 /i _1 Af Q at z ~ 5.7 
(with virial temp erature > lO 4 ^ for gas to cool; see 
iTrac fc Ced 120071) . can also contribute to the two-halo 
term of the stacked profile. According to ITrac fc Cenl 
(2007) (see their Fig.l), the total star formation rate in 
<5 x 1O 9 /i _1 M0 halos can contribute about two thirds 
of the global rate. The amplitude of the clustering with 
these low mass halos is lower, and we expect that in- 
cluding the Lya emission from the <5 x 10 9 /i _1 M & ha- 
los would boost the two-halo term of the surface bright- 
ness profile presented in this paper at the factor of two 
level and that the two characteristic scales would shift 
inward. Our model does not resolve subhalos, but we 
do not expect a large subhalo (or satellite galaxy) frac- 
tion for galaxies at redshift z ~ 6. If we consider a 
mass threshold sample of halos with Mh > M m ; n and 
assume that the mean number of subha los increase lin- 
early with halo mass as M^/(10M m i n ) ()Kravtsov et al.l 
2004), we find that the satellite fraction is about 15% for 
M min ~ 1O u /i _1 M . As subhalos are distributed inside 
virial radius of the parent halo, including star formation 
in subhalos would make the central cuspy surface bright- 
ness profile in the stacked image more extended, which 
would smooth the transition to the plateau part of the 
profile and increase the inner characteristic scale. Inves- 
tigating the exact magnitude of the effect is beyond the 
scope of this paper. 

Besides star formation, there are other sources of Lya 
emission that increase in the dense regions around galax- 
ies, including fluorescence from the UV background and 
cooling radiation. At z ~ 5.7, the fluorescent emis- 
sion caused by the UV background is expected to be 
at the level of a few times 10 ~ 21 ergs~ 1 cm~ 2 a r csec" 2 
(e.g., IGould fc Weinberg! HoM IKollmeier etHI l2f)Tnh . 
much fainter than the signal of interest in this pa- 
per. Cooling radiation may have a noticeable contribu- 
tion to the stacked Lya surface brightness profile, but 
there are uncertainties in the theoretical prediction (e.g., 
iFaucher-Giguere et al.l l2010f) . Further work is needed 
to see whether Lya emission from cooling radiation can 
change the profile and to see how one can separate Lya 
emission caused by star formation from other sources of 
Lya emission. 

Finally, for observational efforts, sky subtraction is an 
important step in revealing the extended Lya emission 
from star-forming galaxies. To compare the observation 
with the theoretical prediction from simulations, a de- 
tailed description of the method of sky subtraction used 
in the observation is necessary, which should be repro- 
duced in the same way in the model. 

4. CONCLUSION 



After escaping from the ISM, Lya photons converted 
from ionizing photons in star-forming galaxies experience 
scatterings in the circumgalactic and intergalactic me- 
dia. Such a radiative transfer process makes Lya emis- 
sion from these galaxies spatially extended. In this pa- 
per, built on the ra diati ve transfer modelin g of LAEs in 
IZheng et all (|2010l) and iZheng et al.l (|201lD . we investi- 
gate the predicted spatial distribution of Lya emission 
that can be measured from the stacked narrowband im- 
age of these galaxies. 

In general, the predicted surface brightness profile 
measured from the stacked image has two characteris- 
tic scales, an inner one at tenths of Mpc and an outer 
one at about 1 Mpc. The profile shows a central cusp 
inside the inner scale, an approximate plateau between 
the two scales, and an extended tail beyond the outer 
scale. The inner scale may possibly be an upper limit, 
as an effect of the grid resolution in the simulation (see 

The stacked surface brightness profile can be under- 
stood as a superposition of the brightness distribution 
from the stacked sources themselves (one-halo term) 
and that from neighboring clustered sources (two-halo 
term). The two- halo term is the profile of the (angu- 
lar) two-point correlation function (plus one) smoothed 
by the extended Lya emission profile (PSF) of individ- 
ual sources. The smoothing makes the profile flattened 
on scales smaller than the spatial extent of Lya emission 
of clustered sources. The outer characteristic scale in 
the stacked surface brightness profile marks this spatial 
extent, and the plateau in the profile is a consequence 
of the smoothing effect. The transition from one-halo 
term domination to two-halo term domination leads to 
the inner characteristic scale seen in the stacked surface 
brightness profile. 

For continuum selected galaxies (LBGs), the amplitude 
of the stacked Lya surface brightness profile increases 
with the source UV luminosity (or halo mass) on both 
small and large scales. The two characteristic scales also 
increase with the UV luminosity. The central cusp be- 
comes steeper for higher UV luminosity. For Lya line 
selected galaxies (LAEs), the amplitude of the central 
cuspy profile increases with observed Lya luminosity and 
the slope is steeper than that of LBGs. Beyond the inner 
characteristic scale, the amplitude of the surface bright- 
ness profile only has a weak dependence on observed Lya 
luminosity. The inner and outer characteristic scales do 
not show strong dependence on Lya luminosity, either. 

Because of the contribution from source clustering to 
the stacked surface brightness profile, the cumulative 
Lya luminosity from the stacked image does not converge 
to the intrinsic Lya luminosity of the stacked sources. It 
is therefore not straightforward to estimate the total Lya 
emission from the stacked image. 

In detail, the surface brightness profile depends on the 
initial line profile of Lya emission (after Lya photons 
escape the ISM). The source would appear to be more 
compact in the narrowband image as the initial line shift 
toward red increases. If this shift is caused by galac- 
tic wind, the wind velocity needs to be comparable to 
a few times the virial velocity of the host halo and the 
wind needs to be largely isotropic to make Lya emission 
point-like. Conversely, from the measured Lya surface 
brightness profile in the narrowband image, it may be 
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possible to constrain the effect of galactic winds. It is 
worth further study along this line. 

Our particular prediction here is for sources after reion- 
ization is complete (z ~ 5.7). It is interesting to investi- 
gate how the prediction changes for star-forming galaxies 
at lower redshifts, because of the large amount of efforts 
in ob serving star-forming galaxies around redshift 2-4 
(e.g., iRauch et al l I2008t ISteidel et all 120111) and in the 
local universe (e.g.. lKunth et al. l l200aiHavesiraIll2007l: 
lOstlin et al J 120 09). It is also necessary to study the ex- 
tended emission for sources at higher redshifts, before 
reionization is com plete, given the increasing observa- 
tiona l efforts (e.g.. lOuchi et all l2010t iKashikaw a et al.l 

At the time of our initial submission, we conclude that 
deep narrowband photometry from large ground-based 
telescopes is on the verge of detecting the extended Lya 
emission aroun d star-forming galaxies, including LBG s 
and LAEs (e.g.. lHavashino et al.ll2004l: lOno et al.H2010D . 
Now the latest observation by ISteidel et al.l ( 201 If ) in- 
deed revealed the extended Lya emission to large radii 
(~ 10") around LBGs and LAEs, starting to verify our 
prediction and provide stringent test to the theoretical 
model. The detection of the predicted Lya emission sup- 
ports the generic picture that Lya radiative transfer in 
the circumgalactic and intergalactic environments pro- 
duces extended Lya emission. The extended Lya emis- 
sion opens a new window to study the circumgalactic and 
intergalactic environment of high-redshift star-forming 
galaxies. The surface brightness profile encodes informa- 
tion about cold baryons around galaxies, including their 



density, temperature, and velocity. All of these proper- 
ties could be modified by galactic wind, so the extended 
emission in principle can put constraints on the galactic 
wind. The stacked Lya image also includes contributions 
from faint galaxies that cannot be detected in a single 
exposure, which provides an opportunity to study low 
luminosity star-forming galaxies. The extended emis- 
sion gives us a better idea on the total amount of Lya 
emission from galaxies. When compared with the UV 
emission or other optically-thin line emission (e.g., Ha 
emission), we may infer the dust distribution and its 
effect on Lya photons. With integral-field-units obser- 
vations of high-redshift star-forming galaxies, we would 
have stacked spectra for the extended emission and ex- 
pect to learn more about galaxy environments. Details 
on how to extract all the information encoded in the ex- 
tended Lya emission need to be investigated, and Lya 
radiative transfer calculation, in combination with so- 
phisticated models of star-forming galaxies and their en- 
vironments, will play an irreplaceable role. 
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